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ABSTRACT 

Radiation from many astrophysical sources, e.g. gamma-ray bursts and active galactic nuclei, is believed to 
arise from relativistically shocked coUisionless plasmas. Such sources often exhibit highly transient spectra 
evolving rapidly, compared with source lifetimes. Radiation emitted from these sources is typically associated 
with non-linear plasma physics, complex field topologies and non-thermal particle distributions. In such cir- 
cumstances a standard synchrotron paradigm may fail to produce accurate conclusions regarding the underlying 
physics. Simulating spectral emission and spectral evolution numerically in various relativistic shock scenarios 
is then the only viable method to determine the detailed physical origin of the emitted spectra. In this Letter we 
present synthetic radiation spectra representing the early stage development of the filamentation (streaming) 
instability of an initially unmagnetized plasma, which is relevant for both coUisionless shock formation and re- 
connection dynamics in relativistic astrophysical outflows, as well as for laboratory astrophysics experiments. 
Results were obtained using a highly efficient in situ diagnostics method, based on detailed particle-in-cell 
modeling of coUisionless plasmas. The synthetic spectra obtained here are compared with those predicted by a 
semi-analytical model for jitter radiation from the filamentation instability, the latter including self-consistent 
generated field topologies and particle distributions obtained from the simulations reported upon here. Spectra 
exhibit dependence on the presence — or absence — of an inert plasma constituent, when comparing baryonic 
plasmas (i.e. containing protons) with pair plasmas. The results also illustrate that considerable care should be 
taken when using lower-dimensional models to obtain information about the astrophysical phenomena gener- 
ating observed spectra. 

Subject headings: acceleration of particles — instabilities — methods: numerical — magnetic fields — plasmas 
— radiation mechanisms; non-thermal 



1. INTRODUCTION 

High-energy astrophysical emission associated with local- 
ized and transient sources - for example gamma-ray bursts 
(GRBs) — is believed to arise from relativistically shocked 
plasmas. The transition from up-stream unshocked plasma 
to down-stream shocked plasma involves a rapid decrease in 
Lorentz factor over the shock ramp, with counter-streaming 
plasmas and associated streaming instabilities as a conse- 
quence - see i.e. Mesz aros] ( |20021 ). Likewise, in outflows 
from relativistic reconnection events, unsteady short lived col- 
lisionless shocks are formed by relativistic streaming instabil- 
ities (i.e. [Zenitani & Hesse] ( |2008 )). Reconnection is often 
associated with flaring events, for example in solar flares and 
in transients of GRBs (flaring/sub-bursts). 

Streaming instabilities, including the filamentation instabil- 
ity (FI), are effectively at work in the propagation of coUi- 
sionless shocks in initially non-magnetized, or weakly mag- 
netized, regions. The instabilities result in the generation 
of complex and time-dependent electro-magnetic fields, with 
structures covering a range of scales. Detailed phase space 
and field topology information is needed in order to ex- 
tract precise diagnostics from such flows, where the stan- 
dard synchrotron paradigm assumptions are explicitly vio- 
lated ^ Medvedev & Loeb |T999)|Silva et al.|2003||Frederiksen| 
let al.| 2004; Spitkovsky|2008|). 

Particle-in-cell (PIC) models are naturally suited for this 
task, since radiation spectra may be synthesized directly from 



such simulations by simply integrating the expression for the 
radiated power, derived from the Lienard-Wiechert potentials 
for a large number of repr esentative particles in the PIC rep- 
resentation of the plasma (|Hededal|2005 Phouc et al. |2005' 
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Nishikawa et al. 2008; Si roni & Spitkovsky|2009p r 

Shock plasma conditions develop through several stages, 
from two-stream non-thermal to almost isotropically ther- 
mal. All radiative stages of the relevant instabilities should be 
probed when modeling optically thin emission spectra from 
highly transient phenomena. Here we exemplify one such 
transition, passing through several radiative stages of the FI, 
with a PIC simulations of a spatially periodic and initially un- 
magnetized two-stream system. 

We emphasize the exploratory nature of this setup, which 
provides a well defined and reproducible model setup where 
such a transition occurs selfconsistently. We make no claim 
that there is a one-to-one mapping of spatial features and time 
scales to real phenomena. 

Section 2 outlines our synthesis method and PIC simula- 
tions while Section 3 is dedicated to results. Discussions and 
conclusions are found in Section 4. 

2. SIMULATIONS 

2.1. Filamentation Instability Simulations 

The relativistic filamentation instability is set up by prepar- 
ing two equal density neutral plasma beams (Fig. [TJ, un- 
der fully periodic boundary conditions. Beams were chosen 
with bulk flows of Tbm e [2,4,6, 10, 15] (Tbm e [4,6, 10]) 
for 3 (2) spatial components and 3 (2) velocity components, 
or 3D3V (2D2V) runs. Simulation volumes were SOOAjc x 
SOOAy X 500Az = 50^6^, (6000A;c x 6000Az = 300^3^), re- 
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Fig. 1 . — Left: transverse magnetic field, Bx(-L.z) along the streaming direction. Yellow arrows show how multiple observers are distributed on Z[obs,z] G 
[0,..,7r/2] w.r.t. the streaming direction, z- Narrow white cones with opening angles S8 = l/Tg,,!;, indicate directions, 9c = Z (z,cohe), 9c S [0, ...,7r/2], of 
solid angle elements chosen in testing for anisotropy in momentum distribution and acceleration of electrons - cf. section [X2] Middle: transverse magnetic 
field, Bx (x, y) across the streaming direction. Two identical plasmas (either + e^ or + p) are counter-streaming, as indicated by the circular symbols. Right: 
temporal development of magnetic field power spectra taken across (B±(kj_)- solid line) and along (Bx(^|| ) - dashed line) the streaming direction z- Black lines 
represents latest simulation snapshot, and colored lines indicate predecessive snapshots. Vertical dashed lines indicate ipeu/t location. 



spectively. Plasma densities were adjusted such that all 
simulations contained an equal number of relativistic elec- 
tron skindepths, 6^ = uj~lc, independent of Thuik, with u 



tegrating the signal as: 



pe 



[A-TTrieq^ /meTbiiikY^^ the relativistic plasma frequency. We 
simulated both baryonic (e" -I- p^) and pair (e~ + e^) plas- 
mas. A reduced mass ratio was used for the baryonic case; 
nii/me = 64 suffices in separating ion and electron dynamic 
time scales. 

In addition, a subset of the radiation spectra were collected 
in static (frozen) EM fields, starting from a well developed 
time in the simulation, rather than in in situ full electromag- 
netic dynamics. These 'frozen' test runs are analyzed and 
discussed by [Medvedev et al.| ( |2010[ ). In summary they con- 
clude that, unlike realistic synthetic spectra collected in situ, 
spectra obtained in static (yet self-consistent) fields are well 
(but erroneously) interpreted by a synchrotron template. This, 
along with the obvious advantage of computational speed and 
ease was a main motivation for employing a completely self- 
consistent in situ synthesis approach. 

We chose the often employed PIC scaling of natural con- 
stants, setting c=mi,=qi,=\. Varying n^, (number or charge 
density) ensured an equal number of skindepths for all TbuII< 
for reasons of comparison. Time is henceforth normalized to 
tw~pl = 1. All results were obtained using a 6''' order field 
solver for Maxwell's equations. The field spectra were only 
affected qualitatively at high wavenumber (i.e. in the extreme 
dissipative range). Spectral synthesis was verified to exhibit 
negligible dependence of field solver order. We used cubic 
(quadratic) spline interpolation between particles and the field 
grid in 3D (2D). 



2.2. Radiation Spectral Synthesis 

We incorporate into the Ph otonPlasma P IC code, devel 
oped in Copenhagen (Haugb0l le |2005[ Hededal 2005), a dis 
cretized version of the standard expression for radiate d power 
from accelerated charges P{uj)(xd^W / dil.duj (e.g. Rybicki & 
Lightman ( |I979 1). Using the aforementioned Pic scaling, 
synthetic radiation frequencies are normalized to . Spectra 
are collected at runtime, and for simulations using A^,o,^10'" 
particles, we collect spectra from Nsy„,i,^lO^ particles with- 
out noticeable slow-down in simulation speed. We retain all 
phase and temporal information of the ensemble signal by in- 
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with 0'=f' - n- ro(f')/c). Spectra are recorded every small 
multiple of a simulation timestep. Thus, we may extract 'in- 
stantaneous' spectral states as P^{ti,Q)- P^{t(),Q). We can 
specify multiple (far-field) observers at arbitrary orientation, 
and we can choose multiple regions for spectral collection in 
the simulation volume. In practice, we have a choice of se- 
lecting phase coherent or phase incoherent spectral synthesis. 
Binning over frequency spectral range was chosen to be log- 
arithmic in Alu, but can in general be either logarithmic or 
linear in Auj; linear frequency binning is used primarily for 
test purposes. Particles selected for synthetic radiation contri- 
bution are subcycled, typically Atsxnth < 10"' Afj,„, to resolve 
high frequency radiation emitted by particles with /3f,„ « I 
and n II /?. Six observers (placed at 'infinity') in 3D (2D) 

are spread uniformly over angles, Z{z, o/js} € [0,7r/2] with 
respect to the streaming direction. Observer "0" designates 
head-on, and observer "5" is edge-on (3D case). 

A number of plasma effects are neglected here. They shall 
be addressed in future work: plasma frequency cutoff, Razin 
effect and self-absorption processes. The spectral synthesis 
method has been sanity checked exhaustively ( Hededal|2005 
.Hededal & Nordlund |2005) . 

3. RESULTS 

3.1. Evolution of the Filamentation Instability 

Growth rates of both linear and non-linear FI depend heav- 
ily on the plasma constituents ( Frederiksen et al.||200 4); in 
e~ + p plasmas, due to the presence of an inert species, 
growth rates are reduced ( [Tzoufras et al.|[2006| l, compared 
with e~ + e'^ plasmas. Emerging electromagnetic field topolo- 
gies and power spectra are likewise affected. For our 
setup two intervals of instability, and a brief turn-over, are 
identified for the relativistic FI: 1) linear, f e [0, 12] (f € 
[0, 14], for e~ + e^ case), 2) saturation/transition, f = [14, 16] 
(f e [12,14]), and 3) non-linear stage, t e [16, "cx)"] (f e 
[14, "oo"]). We shall use these intervals when comparing 
spectra from baryonic and pair plasmas, (electro- )Magnetic 
fields develop rapidly during the linear phase. Subsequent 
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Fig. 2. — panel: total particle distribution (dashed curves) in 7^, above 7^ = Fg^^^j.. Particle distribution subset (solid curves) inside a narrow cone, 
Z{pl,z} < iTg', pointing head-on the streaming direction, i.e. 6c = 0. Right panel: inset "1" plots relative deviation from isotropy of the sub-population- 
to-total population. Curves are offset by factors of X 10 for clarity. Color indicates time; purple (bottom) = earliest, black (top) = latest. Relative deviation is 
defined such that AA'(7^) = 4 = 3 + 1 means three times more particles at given 7^ than is present assuming an isotropic distribution. The 4 narrow panels (insets 
"2"-"5") view the same relative deviation from isotropy now for cones, — F^' < Z{/)t,£} — (m7r/8) < +rg' for m = {1, . . . ,4}, oriented at successively inclined 
angles (indicated by white cones in Fig. [T] leftmost panel) w.r.t. the streaming direction. Color coding designates simulation time, with ; £ [24, 32, . . . , 72] 
(bottom-to-top), respectively, in all panelsT'lots are for the 3D baryonic (e" + p+) case. All plots are Log-Log. NB; This analysis does not assume Maxwellians 
for the (total) isotropic ensemble. We can still have isotropy with non-Maxwellians - which is seen to be the case; even the total ensemble is non-thermal. 



to saturation the field spectra converge to a quasi-stationary 
configuration. A representative case is plotted in the right- 
most panel of Fig. [T] for late times of the evolution. The 
perpendicular Fourier spectra of the magnetic field, Bi^j^ = 
{B ±(k±)B*j_(k±)) , fit a Gaussian -1- a high frequency power- 
law segment on the upper dynamic range, with a dissipative 
cutoff for k^k^y (Nyquist scale). 'Parallel' Fourier spectra of 
the magnetic field, = {B ±(k0B*j_(k^i^)) , are single pow- 
erlaws on < < fc||,wv Except for the earliest linear 
FI, temporal evolution of B± is modest - see Fig. [l] (right 
panel). Particle PDFs are well represented by two shifted 
Gaussians in the streaming direction and a single Gaus- 
sian perpendicular to the streaming direction, f(p7,p±,TB) ^ 
{exp{-(p,-rB)^) + exp{-(p:. + rBf))exp{-p\). The average 
Lorentz factor of the streaming particles decrease in the sim- 
ulation (lab) frame during the FI. Temperatures increase with 
T± > Tji for early times and Tj^wlji in the late phase. Once the 
FI has burned all its momentum anisotropy, the total electron 
ensemble has Thuik ~ 1 (no bulk flow). 

3.2. Non-isotropic electron acceleration during FI 

Non-isotropic momentum distributions naturally result 
from non-isotropic particle acceleration; we expect this to be 
the case based on our devised toy model of "non-Fermi" ac- 
celeration, ( |Hededal et al.|2004[p"edvedev|2006bl l. 

We therefore proceeded to compare subsets of the total par- 
ticle ensemble, in momentum space, specified by the criteria: 

0c,m-S9c < arctan(|/?l|"'/7j|) < dc.m + ^(^c, 

effectively selecting particles with momenta confined inside a 
narrow pitch cone of opening angle S9c = 1 /^Buik, at varying 
inclination, 6'c,m=»J7r/8, m E {0, 1, . . . ,4}, w.r.t. the stream- 
ing direction. Thus, electrons are selected which are signifi- 
cantly beaming their radiative contribution in the direction of 
the cone; they radiate efficiently for certain observer orienta- 
tions, dc- 

Figure |2] (insets 1 & 2) show how, early during the FI, 
electrons experience a non-isotropic acceleration, with direc- 
tional preference towards forward beaming of 1/Tbuik- This 
'beamed' acceleration continues to drive a small part of the 



total particle ensemble to high gamma factors {je.max ~ 45). 
Later, as Tbuik decreases (bi-directionally), this anisotropy 
spreads to include also higher inclination angles (inset 3). 
This could be either due to early 'preference-accelerated' 
electrons diffusing out of their streaming directed pitch cone 
(inset l—> inset 4), or due to current filaments going unstable 
and acquiring a cross-stream component (i.e. they start bend- 
ing), or both. 

3.3. Synthetic radiation, spectral evolution 

Figure [3] reveals consistently higher peak frequencies, 
^peak, for increasing observer angles up to a maximum an- 
gle < \9peak\ < 7r/2. For higher angles ujpeak stays approx- 
imately constant up to 6' ^ 7r/2. Even at late times (Fig. [3] 
lower panel) peak frequencies LUp^ak are observer dependent, 
and apparently increase with streaming axis inclination angle. 
Still, even though LOpeak{\d\ > 0) increases relative to the head- 
on {6 = 0) observer at a given time, it never exceeds the initial 
head-on maximum peak frequency, LUpeakifi = 0), in 3D. 

In the 2D case viewed in Fig.|4]the peak frequency increases 
for all observers during the FI, and cUpeak.iD > ^peak,3D for all 
cases (open circles of both Fig.|3]and Fig.|4]i. This effect is an 
artifact of the lower-dimensional 2D model; particles experi- 
ence more frequent interactions with the ion filaments, since 
the "filament hit probability" is considerably smaller in 3D 
due to the reduction in filament filling factor (filaments are ac- 
tually 'sheets' in 2D, whereas they are true line-like filaments 
in 3D). Radiation spectra are therefore formed more rapidly 
by the plasma field structures and with better statistics in 2D 
- albeit artifactually so. 

Fitting the semi-analytical result for jitter radiation 
(Medv eHev et al.|2010( l to 3D spectra (Figure |5] upper panel) 
shows markedly better fits than a synchrotron spectrum for 
the same parameters, during early stage FI, both for bary- 
onic plasmas and pair plasmas. At later stages, also mod- 
els of synchrotron fit the spectra well, but we emphasize pre- 
vious findings, that the underlying plasma dynamics is not 
synchrotron in origin Medvedev (2000, 2006a). For a spec- 
trum (and evolution) which has been isotropically averaged 
over all observer directions — Figure |5] lower panel — mak- 
ing the same qualitative scaled fits of the semi-analytic model 
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produces better agreement if a jitter radiation environment is 
assumed. Such an averaged spectrum can be seen as the im- 
mediate spectral state in the bulk of the plasma two-stream (or 
locally in a shock), when the simulation frame coincides with 
the co-moving frame of the plasma, as is the case in the FI 
setup presented here. 

Our result farther agrees quite precisely with the scaling 
relation of peak frequencies for a semi-analytical jitter radi- 
ation model, reported in ( [Medvedev e t al. 2010). At early 
times, we have measured ujpeak.iD ^luik^ while at late times 
ijJpeak,3D oc T^^^j. joint by a smooth transition. 

A 'by-eye' re-alignment fit (thick grey line. Fig. |4|i of the 
same semi-analytic result to the 2D baryonic case shows an 
under-estimate of high frequency radiated power. This under- 
shoot artifact is not seen in 3D (upper green, solid, curve. Fig. 

It, too, is probably a consequence of the reduced dimen- 
sionality, again by the above argument concerning 2D result 
artifacts. 

Moreover, in 20, the maximum value corresponding to the 
abscissae of Fig. [2] -fl^^^^ w 2100, is jf^ 2D ~ 4200 which 
is higher by a factor of 1 max. id hi max 3D 2. This further influ- 
ences late time spectra (lower panels of Fig.[3]and Fig.|5]l. The 
2D and 3D cases exhibit gross discrepancy in peak frequen- 
cies; they differ by an order of magnitude (uj peak, id peak .3D ~ 
10), due to a general upward drift of ujpeak in 2D. Again, 2D 
spectra are pushed to higher frequencies, likely due to an en- 
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Fig. 3. — 3D: Spectra for multiple observers. Upper panel: Early stage 
FI. Spectra have barely had time to form. Lower panel: Same but for late 
stage FI. Open circles (red) added, suggesting variation of peak frequency 
with observer angle. Dashed and dotted lines suggest low and high frequency 
powerlaw slope indices, a ("a") and /3 ("b"), respectively. /3 indicates range 
of possible slopes. 
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Fig. 4. — 2D: Spectra for multiple observers. Upper panel: early stage FI. 
Open circles (red) and thin broken line (black) added to suggest the variation 
of iiipf„j^ with observer angle. The off-axis variation of it^pgak is more pro- 
nounced in 2D; the probability that particles do not pass through a filament 
is much smaller than for 3D. Additionally is given a "by-eye" fit (grey thick 
line), by realignment of the semi-analytical model (cf. the true 3D fit Fig.|5}. 
Dashed and dotted lines have been added to suggest low and high frequency 
powerlaw slope indices, ol ("a") and /3 ("b"), respectively. Lower panel: Late 
stage FI. High frequency slopes are consistently /3 ~ — 1.5 which is harder 
than 3D for all observers except head-on (cf. Fig. (3) . This is consistent with 
a 'filament hit probability' argument (cf. sect. |3 j| 

hanced 'filament hit probability'. The 2D results for + e~- 
plasmas exhibit analogous features. 

4. DISCUSSION AND CONCLUSION 

The context of the current investigation is to compare the 
jitter and synchrotron radiation approaches to radiation spec- 
tra modeling in a simple but well defined and reproducible 
setup. The more general problem of synthesizing the radi- 
ation field from fully developed collisionless shocks, with a 
firm tuning of parameters, is the subject of future studies. 

Here we have produced synthesized radiation spectra for 
multiple observer orientations from 2D and 3D particle-in- 
cell simulations of the relativistic filamentation instability in a 
periodic box, relevant to relativistic collisionless shocks such 
as believed at work in for example gamma-ray bursts, and in 
transient shocks at reconnection sites. We investigated radi- 
ation spectral evolution for both pair plasmas (e" -I- e~) and 
baryon loaded (e~ +p) plasmas, for a broad range of bulk flow 
speeds. 

In summary the results are: 

• At early times, LUpeak,3D oc T^, ,;,, which fits a semi 



analytical j itter radi ation model (Medvedevetal. 20101. 
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Fig. 5. — 3D, baryonic plasma, head-on observer. Upper panel: compar- 
ison of synthesized spectra with semi-analytical model fits to jitter (green, 
solid) and synchrotron models (red, dashed). Spectra for t S [6, 16] and 
t G [18,28], for linear and late stages of the FI, respectively. Late stage is 
offset by factor of xO.2 for clarity. NB; the semi-analytical fit in the bary- 
onic case is "by-eye fit" of that for the pair plasma case, but realigned. The 
functional form still provides the better fit of the two models (jitter vs. syn- 
chrotron). Inset: same fits as in baryonic case, but for a pair plasma, with 
late stage spectrum offset by X 0. 1 . This strict theoretical fit aligns remark- 
ably well with the jitter case. Lower panel: same as upper panel except for 
an isotropically averaged spectrum based on 6 observers at inclination angles 
w.r.t. the streaming direction as described in section [2!2] 

Later ujpeak,3D oc ^luik' with a smooth transition between 
the two regimes. We emphasize that this scaHng is not 
synchrotron in origin. 

It may be difficuh to extract such scaling-laws di- 
rectly from existing observations. However, the result 
is still quite relevant: through a shock interface, the 
bulk Lorentz factor and spectral parameters vary dras- 
tically. In employing semi-analytical models based on 
synchrotron modeling, inference of a most likely physi- 
cal scenario might fall short, yielding incorrect scalings 
of frequencies and flux. 

• Non-isotropic acceleration occurs during both linear 
and early non-linear FI, and is present in both bary- 
onic and pair plasmas. For pairs the process pro- 
gresses much faster, and non-isotropic acceleration is 
negligible. The findings are commensurate with pre- 
vious findings on direct 'non-Fermi' particle acceler- 
ation (Hededal et al. 2004) in non-hnearly developed 
e~+p (i.e. inert) filaments ( Frederiksen et al.|20041 l. 

• Synthetic spectra generated from 2D PIC simulations 
exhibit an upward drift of ujpeak.iD which is not seen in 



• Further, generally, 2D spectra differ qualitatively and 
quantitatively from 3D spectra in both ujpeak, Ppeak = 
Pip^peak) and in spectral slopes, a {(3) at low (high) fre- 
quency - see Fig. [3] and Fig. [5] 

We find that lower-dimensional models (such as for exam- 
ple 2D2V or 2D3V) provide distorted measurements of pa- 
rameters central to radiation modeling. Peak frequencies are 
generally too high, and low- (high-) frequency spectral slopes, 
a (/3), are too steep (flat) compared with spectra obtained 
from 3D modeling. We expect that 2D3V (2D spatial dimen- 
sions, 3D particle velocity components) simulations will dis- 
play properties intermediate of 2D2V and 3D3V. 

While obtaining synthetic spectra from both pair plasmas 
(e^ -I- e") and baryonic plasmas (e" -I- p), we establish that 
spectra are qualitatively dependent on the plasma constituent 
species. This difference must be addressed when interpret- 
ing astrophysical observed spectra. We interpret the radiation 
spectra differences between baryon and pair plasma cases as a 
direct signature of screened current channel formation in the 
non-linear stage of the FI. 

The presence of non-isotropically distributed electrons, ac- 
celerated preferentially in the streaming direction, could lead 
to spectra that peak at higher frequency than would be esti- 
mated from the isotropic (thermal -1- powerlaw) particle en- 
semble. Since they are powerlaw distributed (left panel. 
Fig. [2]) their contribution should also be powerlaw by origin. 
Although non-isotropic acceleration is also present in e^e~- 
plasmas it is negligible on long timescales, due to the lack of 
an inert plasma species. We conjecture that baryonic plasmas 
could therefore yield radiation signatures with harder power- 
law behavior at high-frequency ifibaryonic > Ppair), but softer 
spectral slope at low frequency (aharyonw < cupair)- The elec- 
trostatic field, although screened, also plays a significant role 
for spectrum of emitted radiation. In the up-stream region 
ahead of the shock, where streaming instabilities are of early 
stage type, screened electrostatic effects are likely to produce 
heavily observer dependent effects. For global shock sim- 
ulations with non-periodic boundary conditions such a non- 
isotropic acceleration mechanism will be sustained. 

Our results establish tractability of an alternative approach 
to radiation modeling - the jitter radiation approach. It con- 
trasts the standard synchrotron approach; assuming a pure 
synchrotron mechanism in observed relativistic streaming 
phenomena, we would infer — incorrectly — that large scale 
homogenous magnetic fields were present, when rather it 
might be small scale fields, in (a priori) non-magnetized rela- 
tivistic shocked outflows. 

A tight correspondence exists between on one hand the 
spectral evolution of radiation emitted from electrons in op- 
tically thin coUisionless plasmas and, on the other, the field 
spectral evolution in the plasma due to the instabilities that 
create and mediate such radiation. 

Care should be exhibited with reduced dimensionality mod- 
eling; as 3D models become increasingly affordable they are 
to be preferred over 2D models, as spectral diagnostics can 
be distorted as a consequence of reduced dimensionality — 
filament entanglement is for example only possible in 3D. At- 
tention should also be given to choices of plasma constituents 
(electrons, protons, ions, etc.). 

To conclude, the radiation spectral modeling results pre- 
sented here reveal limitations that need to be considered while 
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reverse engineering spectra observed from, for example, rel- 
ativistic astrophysical outflows, GRBs, relativistic reconnec- 
tion sites, and from laser-plasma interaction in laboratory as- 
trophysics, when examining likely physical scenarios for their 
match to observed spectra. 

AN, JTF, and TH acknowledge support from the Dan- 
ish Natural Science Research Council. MVMs work has 



been supported by NSF grant AST-0708213, NASA ATFP 
grant NNX-08AL39G and DOE grant DE-FG02-07ER54940. 
MVM also acknowledges support from The Ambrose Monell 
Foundation (IAS) and The lb Henriksen Foundation (NBIA). 
Computer time was provided by the Danish Center for Sci- 
entific Computing (DCSC). This work was supported by the 
European Commission through the SOLAIRE Network (AN). 



REFERENCES 



Frederiksen, J. T., Hededal, C. B., Haugb0lle, T., & Nordlund, A. 2004, ApJ, 
608, L13 

Haugb0lle, T. 2005, PhD thesis, Niels Bohr Institute |astro-ph/05 10292 1 
Hededal, C. B., Haugb0lle, T., Frederiksen, J. T., & Nordlund, A. 2004, ApJ, 
617,L107 

Hededal, C. 2005, PhD thesis, Niels Bohi- Institute | astro-ph/0506559J 
Busk Hededal, C, & Nordlund, A. 2005, arXiv:astro-ph/ 05ff662| 
Medvedev, M. V.; Loeb, A., ApJ, vol.526, 697-706 (1999) 
Medvedev, M. V. 2000, ApJ, 540, 704 
Medvedev, M. V. 2006, ApJ, 637, 869 
Medvedev, M. V. 2006, ApJ, 651, L9 

Medvedev, M. V., Frederiksen, J. T., Haugb0lle, T. & Nordlund A, 2010, 

submitted to ApJ, preprint at arXiv: 1003.0063y 2 [astro-ph.HE] 
Meszaros, P 2002, ARA&A, 40, 137 

Nishikawa, K. I., et al. 2008, conf. proc.,"Blazar Variability across the 
Electromagnetic Spectrum", Palaiseau, France (2008) 



Ta Phuoc, K., Burgy, F, Rousseau, J.-P , Malka, V., Rousse, A., Shah, R., 
Umstadter, D., Pukhov, A., & Kiselev, S. 2005, Phys. of Plasmas 12, 
023101 

Rybicki, G. B. & Lightman, A. P., "Radiative processes in astrophysics", 

Wiley & Sons, NY, 1979 
Silva, L. O., Fonseca, R. A., Tonge, J. W., Dawson, J. M., Mori, W. B., & 

Medvedev, M. V. 2003, ApJ, 596, L121 
Sironi, L., & Spitkovsky, A. 2009, ApJ, 707, L92 
Spitkovsky, A. 2008, ApJ, 673, L39 

M. Tzoufras, C. Ren, R S. Tsung, J.W. Tonge, W. B. Mori, M. Fiore, R. A. 

Fonseca, L. O. Silva 2006, Phys. Rev. Lett., 96, 105002 
Zenitani, S., Hesse, M. 2008, Phys. of Plasmas 15, 022101 



